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NOMENCLATURE 


Hertzian contact radius, m 

Specific heat J/(kg K) 

-3 2 

a constant = 10 Ws/m 

Characteristic fluid temperature, C 

Film thickness, m 

Thermal conductivity, W/ (inK) 

2 

Pressure, N/m 

Temperature-viscosity coefficient 

Temperature, C 

Traction coefficient 

Thermal diffusion time, s 

Characteristic fluid resident time, s 

Velocities in the directions, ra/s 

Sliding velocity, m/s 

Normal load, N 

Cartesian coordinates, m 

2 

Viscosity, Ns/m 

Ratio of film thickness to composite surface roughness, h/cf 

3 

Density, kg/rn 

/2 2 

Composite surface roughness = 

2 

Shear stress, N/m 

Dummy integration variable 



subscripts 


1 Stationary surface 

2 Moving surface 

b Ball surface 

c Contact center value 

f Fluid 

H Hertzian conditions 

o Ambient conditions or minimiun value 

s Sapphire surface 
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I . SUMMARY 

During the past year research in the areas listed below has progressed 
and is reported herein: 

Prediction of Traction in Sliding EHD Contacts 
The shear stress theory proposed in the previous report when applied 
to the sliding EHD point contact using a naphthenic oil was found to yield 
results consistent with experimental findings for the high end of the speed 
range studied. The particular speed at which this theory yields reasonable 
results depends on the normal load and film thickness. At lower speeds ^ 
the calculated tractions exceeded the measured values. This is possibly due 
to either the onset of asperity interaction as the film thickness decreased 
or the breakdown of the Newtonian fluid model at high shear stress. 

Calculations based on the theory indicated that the traction tends 
to increase with al^ increase in base viscosity and a decrease in temperature 
viscosity dependence of the lubricant. A weaker traction increase resulted 
from an increase in the pressure-viscosity dependence of the fluid. Varia- 
tion in the thermal conductivity was found to have very little effect on 
the traction . 

Experimental Investigations in Sliding EHD Contacts 
A complete mapping of the average fluid and ball surface temperature 

2 

has. been obtained for the relatively high Herts pressure level of 1.51 GN/m . 

. 2 
The previous ■ report contained similar data for 1.05 GK/m: Hertz pressure. 

These data are compared to yield a qualitative dependence of fluid tempera- 
ture levels on Hertz pressure. Further studies of the bail surface 

2 

temperature only were made at pressure levels from 0^521 to 2.03 GN/m . 



The results of these studies are then compared with the flash temperature 
theories of Blok, Jaeger and Archard, Deviations from the predictions are 
then discussed. 

A detailed study of the effect of surface roughness on EHD contact 
behavior has also been undertaken. Balls of three roughness values 
ranging from .011 to .381 ym c.l.a. were studied using both a surface 
profilimeter and the EHD contact simulator. Ball surface profiles were 
determined prior to running in the EHD contact, and then through the 
use of a relocation technique, the surface profiles were again measured 
after an EHD experiment. It was found that changes in the profiles, as 
a result of asperity interaction, were closely related to the ratio of 
the film thickness to composite surface roughness and the surface tempera- 
ture level. 

Finally, a technique has been developed which allows the determination 
of temperature fluctuations arising from individual asperity interactions. 
All previous temperature data has been time-averaged by the infrared 
detector. The high frequency fluctuations are shown to exist only in 
cases where asperity interaction occurs. It has also been found that the 
range of these fluctuations increases v?ith a decreasing ratio of film 
thickness to composite surface roughness. 
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II. THE PREDICTION OF TRftCTION IN SLIDING EHD CONTACTS 
A. Introduction 

The determination of the traction force in elastohydrodynamic 
contacts is of primary importance in the understanding of many lubricated 
mechanisms. The traction is defined as the force generated in the contact 
which resists relative motion of the bearing surfaces. This is directly 
related to the power loss in mechanical components such as gears, cams, 
and rolling element bearings. Traction studies are also motivated by a 
need to predict the onset of skidding in rolling element bearings. How- 
ever, due to the physical complexity of the problem, no simple model 
describing quantitatively the generation of traction in an EHD contact has 
yet been generally accepted- Complications caused by the large role of 
thermal effects and by difficulties in finding a rheological model which 
is adequate at high pressures and shear stress have hindered simplified 
theoretical analysis - 

A number of models for predicting traction have been proposed. 

Crook Cl) , Kannel and Walowit (2), and Allen, Townsend and Zaretsky (3) 
all assumed isot .ermal. bearing surfaces in their models. Based bn tempera- 
ture measurements by Turchina, Sanborn and Winer {A) this condition appears 
to be physically violated, particularly when a great deal of sliding is 
present. Cheng and SternTicht C5) and Cheng (6) included thermal effects 
in a numerical analysis dealing primarily with line contacts in the pre- 
diction of film thickness, pressure, and temperature. ?^le these inves- 
tigations give a great deal of insight into the effects of thermal behavior 
in the analysis, the complexity of the numerical iteration technique 
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limits the practical application of this model for traction prediction 
over a wide range of physical situations. Trachman and Cheng (7) predicted 
traction in line contacts in a numerical solution which included hoth 
thermal effects, and non-Newtonian rheological models, hrchard and Baglin 
(8) assumed isothermal conditions, a Newtonian fluid, and an exponential 
pressure-viscosity behavior. Because of these assumptions, their model is 
limited to low sliding speeds where thermal effects are less important. 

Ar chard and Baglin derived non-dimensional groups from physical considera- 
tions and established a foundation on which solutions with thermal effects 
and different rheological models may be built, 

Jakobsen and Winer (9-11) have taken thermal effects into account, 
while limiting to a few graphical steps the work required to find the shear 
stress at a point in a sliding EHD contact. The method is thus suitable 
for predicting tractions over a wide range of conditions while eliminating 
the restrictions imposed by the assumption of isothermal walls. In addi- 
tion, Jakobsen and Winer (9,10,12) have proposed a Newtonian model for 
lubricant behavior, employing a power-exponential temperature-viscosity 
function. This model resulted from the correlation of measurements on a 
capillary viscometer at shear stress levels only 3-5 times less than the 
average shear expected in an elastohydrodynamic contact. 

This study attempts to evaluate the shear stress theory and liabri- 
cant rheological model of Jakobsen and Winer C9-12) by applying them to 
traction calculation in a sliding elastohydrodynamic point contact, and 
comparing the caicul at ed values to experimental measurements. The correlation, 
is limited to a single naphthenic oil N1 (15). The theory is also used 
to predict the effect bn traction of variations in the lubricant Tnaterlal 
properties. 
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B. A pplication of the Theory 

The assumptions which form the basis for this theory are summarized 
here in the context of their application to a sliding EHD contact. The 
flow geometry, with major assumptions, is shown in Figure 1. surface 1 is 
stationary, and surface 2 moves with constant velocity in the direction. 
According to both theory and experiment, the change in film thickness over 
most of an EHD contact is very small compared with the size of the contact. 
The bearing surfaces are therefore assumed parallel in the vicinity of the 
point at which the shear stress is being calculated, although the film 
thickness is allowed to vary from point to point. P?hile temperature vari- 
ations are recognized to exist in the x^ and x^ directions, the gradients 
of temperature in these directions are smaller than in the direction by 
about four orders of magnitude. Therefore, temperature is assxamed to be 
pointwise independent of x^ and x^. 

Because surface 2 is moving, the major portion of heat generated 
in the contact is carried away by this surface v^hich justifies the assumption 
that the stationary surface is adiabatic. Over the major portion of the 
contact, the component of shear stress due to the pressure gradient is 
much smaller than that due to sliding motion. Consequently, the pressure 
is assumed pointwise constant. Detailed discussion of these assumptions 
is found in Jakobsen C9) . In addition, body forces are neglected, and 
density and thermal conductivity of the lubricant are assumed constant. 

A Newtonian rheological model of the lubricant, employing a linear 
relationship between shear stress and shear rate, is used. The yiscosity- 
temperature dependence is expressed by the power-exponential relation 
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Assumptions 

K = u^{x2); = V 

2. = U 3 =0 

3. T = Tl*,)' T(0) » T ; T(h) = 

4. constant pressure 

5. constant density 

6. constant thermal conductivity 

7. steady state 

8. laminar flow 

9 . body forces negligible 

Figure 1. The Flow Geometry and Major Assumptions. 
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n = c exp I (s'/T) ^3 (1) 

both E and Q are functions of pressure. Jakobsen (9) has shown that this 
equation gives a reasonable description of lubricant viscosity within the 
range of temperature and shear stress normally found in an elastohydrodyn" 
aiTiic contact. 

Through integrations of the energy equation and the equations of 
motion, reduced by the above assumptions and the given rheological model/ 
the theory provides a means of calculating the shear stress T at any point 
in a sliding EHD contact. The condition of an adiabatic stationary wall 
yields the relation 

f ^2 -1 

ri{T,Q) dT = -^ (2) 

which is used to determine the temperature of the stationary surface, T^. 
The shear stress at a point in the contact is then determined from the 
equation 



f rTicg,Q)"^dg3”^/^ dT 

T 

1 


C3) 


The following parameters must be known: the temperature of the moving 

surface, T^; the film thickness, h; the thermal conductivity of the fluid, 
k; the sliding speed, y; and the visoosity parameters Q and E. These 
latter taro are calculated from the above equation and from viscometric 
data on the lubricant at the pressure of the point under, consideration. 
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In order to perform this calculation, therefore, some assumption about 
the pressure distribution in the contact must be made. A Hertzian 
pressure distribution was chosen because of its simplicity for practical 
application , and its apparent applicability to a wide range of operating 
conditions. The computer^generated solutions to the BHD problem of Cheng 
and Sternlicht (5) show pressure distributions which converge to the 
Hertzian as the speed decreases and the load increases. 

The procedure for calculating the shear stress at a point in the 
EHD contact is represented in the form of dimensionless charts in Jakobsen 
and Winer (11). The basic equations of the present theory were solved 
numerically. The traction in an EHD point contact was found by systemati- 
cally calculating the shear stress throughout the Hertzian contact, and 
numerically integrating it over the contact area to yield the traction. 
Details of the numerical methods used, as well as a listing of the computer 
program, are given in Kuns (13) . 

C. Evaluation of the Theory 

Measurements of traction for mineral oil N1 were made in a sliding 
EHD point contact using the equipment described in Sanborn and Winer (14) » 
and the results compared to calculations using the methods described in 
the previous section. The conta-t consisted of a rotating steel ball 
whose surface was sliding with constant speed on a stationary flat surface 
of synthetic sapphire. The traction was measured for normal loads of 67 
and 334 H, and sliding speeds ranging from .35 to 12.7 m/s. . 

The temperature distribution on the- surface of the moving ball 
for each set of conditions was known from measurements made using the 
infrared technique described in Turchina, Sanborn and Winer (4). The 
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distribution of film thickness in the contact area was also known from 

measurements made using the optical interference method of Sanborn and 

Winer (15) . The required pressure-viscosity data for the fluid N1 was 

obtained from Novak (IS) . Because viscosity data vras only available at 

8 2 

pressures up to 1.38 x 10 N/m ^ it was necessary to extrapolate the 

9 2 

experimental curves up to the maximum Hertzian pressure of 1.0 x 10 N/ra 

9 2 

for the 67N load, and 1. ’ x 10 N/m for the 334W load. This extrapolation 
was performed using Roelands’ (17) pressure-viscosity correlation. 

Figure 2 shows a plot of both calculated and measured tractions 
against sliding speed for the 67 N load. The traction is represented in 
the form of a traction coefficient, defined as the ratio between the trac- 
tion force and the normal load. Figure 3 is a similar plot for the 334 N 
load. Figures 2 and 3 both indicate that the agreement between calculated 
and measured traction becomes better as the sliding speed increases. 

Because the film thickness decreases as the sliding speed decreases, there 
is a possibility that the divergence at lower speeds is related to asperity 
interaction. Such interaction could affect the contact in basically two 
ways which would tend to decrease the actual traction below the calculated 
values, as is observed in the figures. First, the contact between asperi- 
ties could act as an additional source of heat which would raise the film 
temperature, and in turn lower the lubricant viscosity and therefore the: 
traction. The traction model used in the calculations assumes that the 
only source of heat in the EHD contact is that due. to viscous shearing 
of the lubricant. Therefore, one could expect calculated tractions to 
exceed the measured values. In addition, as the film thickriess decreases 
and asperities come into contact, the asperities themselves would be 
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carrying a part of the load. This would result in lower pressure on the 
fluid and therefore lower viscosity and a lower measured traction than 
predicted. Asperity interaction could cause traction to increase with 
some surfaces, hut in the present case the surfaces were hard (52100 
against sapphire) which minimized plowing and shearing effects. 

In order to further investigate the possdhility of asperity inter- 
action, surface profiles of the steel ball and the sapphire were recorded 
on a profilometer. A composite rms surface roughness of about .046 yra 
was found for the two surfaces. A study by Sibley (18) indicates that 
some asperity interaction may be expected when the ratio A. of film thick- 
ness to composite surface roughness is less than approximately 2.5. Recent 
surface temperature measurements in this laboratory (see Section III) 
indicate that asperity interaction begins to influence surface temperatures 
as A becomes one or less. Accordingly, the measured film thickness at 
each speed is shown in Figures 2 and 3, along with the region cross hatched 
which indicates a A of from 1 to 2.5. if the film thickness is below the 
cross hatched region, asperity interaction is certain to occur. 

. Another consideration in explaining the disparity between calculated 
and measured results is the behavior of the fluid under the extreme condi- 
tions imposed by a low sliding speed. It may be unreasonable to assume that 
a single rheological moi3el is adequate to describe lubricant behavior 
throughout such a wide range of operating conditions. Johnson and Roberts 
(19) suggest that, above some critical shear stress, the fluid film behaves 
in the manner of a plastic solid rather than a viscous liquid. The possi- 
bility therefore exists that the rheological model used in this work 
breaks down when the shear stress in the film reaches some maximum 
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limiting value. The Johnson and Roberts limiting shear stress model 

would result in tractions closer to the measured values than the present 

model if the limiting shear stress of the fluid is less than the values 

calculated with the present model. 

For the 67 W load, the maximum shear stress calculated at 1.40 ra/s 

7 2 

was in the neighborhood of 8.0 x 10 ^l/m , Similarly, for the 334 N load at 

2.54 m/s, the maximum shear stress in the contact was calculated to be 
7 2 

8.1 X 10 W/ra . Each of these two cases corresponds to the transition 
speed at which calculated and measured tractions began to deviate as the 
speed decreased, at higher speeds/ calculated shear stresses were lower 
than these values throughout the contact- At speeds lower than the 
transition speeds, these shear stresses were exceeded by the calculated 
value. The fact that these maximum shear stresses are essentially the 
SEane for the two loads lends credence to the possibility that the Nev/tonian 
rheological model of the lubricant used in the calculations fails to 
adequately describe fluid behavior at extremely high shear stresses. 

D. The Effect of Material Parameters on Traction 

The effects on the traction of variations in the lubricant material 
parameters; temperature-viscosity dependence, pressure-viscosity depen- 
dence, base or inlet viscosity, and thermal conductivity were studied. 

These quMti ties were Varied independently in calculations using the 
measured surface temperature cind f lira thickness found for Fluid N1 {the 
naphthenic base oil {15} 1 . The results : were then compared with the cal- 
culations for Fluid Wl. 

3ji order to examine the effect of teraperature-viscosity dependence , ■ 
on the calculated traction, a model fluid having the same pressure dependence 
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and inlet viscosity as Fluid Wl, but having an exponential temperature- 
viscosity coefficient I— one-half that of N1 at inlet conditions 

was used in the program. The viscosity-temperature behavior of this model 
fluid is similar to that of a dimethyl siloxane, a lubricant with low 
temperature-viscosity dependences but the pressure-viscosity chfiracteristics 
are quite different. The resulting calculated traction coefficient for a 
67 N load and a speed of 2.54 m/s was three times greater than that cal- 
culated for Fluid ISJl for the same load and speed conditions. The trac- 
tion is therefore increased by a decrease in the temperature-viscosity 
coefficient of the lubricant, and appears to be fairly sensitive to such 
changes. 

The pressure viscosity dependence was investigated by using a 

different means of extrapolating Novak's (16) viscosity data for Fluid Nl 

to the high pressures near the center of the contact. The shape of the 

pressure-viscosity curve was taken to be the same as that for a similar 

naphthenic fluid, (Fluid 36-G of the ASB-IR Pressure-Viscosity Report (20) 

3 2 

for which data was available up to 6.9 x 10 N/m . The viscosity thus 
obtained for the pressure at the center of the contact was three times 
less than that previously used for Fluid Nl. The viscosity and temperature 
dependence at inlet conditions were taken to be the same as those for 
Fluid Nl « The traction calculations using these viscosities were about 
15% lower than those of Fluid Nl for the same load and speed conditions. 
This indicates that an increase in the traction would result from increases 
in the pressure-viscosity dependence,, but that the traction is relatively 
insensitive to such changes. 

The effect of the inlet viscosity on the traction iras also studied. 
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another model fluid was used, having the same pressure and temperature 
dependence as Fluid Nl, but having a viscosity at inlet conditions three 
times that of Fluid Nl. as before, all other system parameters were 
held constant. The resulting traction coefficient was greater than that 
for Fluid Nl by a factor of 1.75. The traction therefore increases as 
the inlet viscosity increases. 

The only other lubricant material parameter involved in the traction 
calculation is the thermal conductivity. Variations in this property appear 
to have very little effect on the traction. In addition, the thermal 
conductivities for most hydrocarbon and dimethyl siloxane lubricants fall 
within a very narrow range. 

TT'action measurements would not be expected to be as sensitive to 
material property differences as the above calculations indicate. This is 
primarily because the surface temperatures for Fluid Nl were used in the 
calculations. Using a more viscous fluid, for instance, would increase 
the energy dissipation in the contact, thus increasing the temperature 
and decreasing the traction until a lower steady state value is reached 
with higher surface temperatures. Temperatures and film thicknesses for 
Fluid Nl were used only so as to allow independent variation of the 
material parameters. Consequently, the results of the calculations of this 
section should be viewed as relative variations only, rather than as 
absolute magnitudes to be expected in actual lubrication situations. 

E. Conclusions 

The shear stress theory of Jafcobsen and TSiner (11) , when applied 
to a sliding EHD point contact, was found to yield realistic values for 
the traction at high sliding speeds. The particular .speed at which the - ^ 
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theory becomes applicable depends on such factors as normal load and film 
thickness, at lower sliding speeds , the calculated tractions exceeded 
the measured values. This was possibly due either to the onset of asperity 
interaction as the film thickness decreased, or to the breakdown of the 
Newtonian fluid model at high shear stress. 

Calculations based on the theory indicated that the traction tends 
to increase with an increase in base viscosity and a decrease in tempera- 
ture viscosity dependence of the lubricant. A weaker traction increase 
resulted from an increase in the pressure-viscosity dependence of the 
fluid. Variation in the thermal conductivity was found to have very 


little effect on the traction. 
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III, EXPERIMENTAL INVESTIGATIONS IN SLIDING EHD CONTACTS 

A. Temperature Distributions at a 215N Normal Load 

The distribution of the average fluid film temperature and ball 

2 

surface temperature, for a load of 67 N (1.05 GN/m maximum Herts pressure) 

have been obtained using the test apparatus shown in Figure 4 and the 

results were reported previously [10,21]. Similar measurements for the 

2 

case of a normal load of 215N (1.51 GN/m maximum Herts pressure) were 
made and the most important results are given in Table I. It can be seen 
from Table I that the film thickness at the center of the .EHD conjunction, 
at the side lobes and at the contact exit are all essentially the same. 
This constant film thickness tends to give more symmetric temperature 
distributions than previously observed. 

In addition, at the higher load (215N) , the maximum ball surface and 
fluid temperatures always occurred on or near the contact centerline at 
a point downstream of the contact center. The reason for this is that 
the pressure, and thus viscosity, is greatest at the contact center, and 
consequently,, for a uniform film thickness, viscous dissipation is a max- 
imum at the center. However, since the fluid residence period increases . 
as it moves toward the exit, the position of the maximum temperature 
could be between the contact center and exit depending on the effective-? 
ness of heat transfer at the bearing surfaces. At lower loads [10] , the 
maximum contact temperature may occur off 1die contact center line, for 
example, at the side lobe constriction. This is attributable to the very 
low film thicknesses at the side lobes compared to other locations. 

Also, the position of the maximum fluid film temperature is upstream 
of the position of the maximum ball surface temperature. Keeping in mind 



Table I. Suimnairy of Experimental Results for 215N Normal Load 
(1.51 GN/m^ = 219 kpsi peak Pressure) 
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(Eluid Nl, Load 215N (peak Pressure = 1.51 GN/m^) , Bath Temp. 40+ 1C). 
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Figure 4. Sliding EHD Test Apparatus. 
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that the fluid film temperatures are essentially fourth power averages 
through the film, the observed shift in peak temperatures does not imply 
a violation of physical principles since the fluid temperature gradients 
at the ball surface may change considerably with position along the center- 
line allowing a matching of ball and fluid temperature levels at the 
interface. 

Figures 5 and 6 show temperature contours for a sliding speed of 
1.39 m/s. The general shape of the ball surface temperature contours 
(Figure 5) is similar to that found at lower loads [10] . In the present 
case, the ball surface enters the Hertaian contact at 85 C and reaches a 
maximum of 180 C downstream of the contact center. The fluid temperature 
contours shown in Figure 6 show that the difference between the primary 
maximum (226 C) on the centerline and the secondary maximum (186 C) in 
the side lobes has significantly increased compared to the lighter load 
of 67N [10] . In the lighter load case, the temperatxires were 171 C and 
173 C respectively. At 215N, the secondary maximum disappears as the 
speed is increased. 

Figures 7 and 8 show respectively the ball surface and fluid tempera- 
tures along the centerline for a range of speeds. Results for the 67N 
load are also shown. The difference between the nature of the temperature 
distribution along the centerline for the two loads is attributable to 
different shapes of the film thickness profiles. The nearly constant film 
thickness at the 215K load yields fluid temperatures, which like the ball 
temperatures in both cases, change smoothly along the centerline ^ 

Figure 9 shoves the inlet temperature increase as. functions of speed 
for the fluid and ball surface at the higher load condition. Temperatures 
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have been measured at the point of entry to the Hertaian contact zone. 

It is apparent that a significant amount of inlet heating has occurred. 
Inlet shear heating is important because it has a significant effect on 
the viscosity of the oil entering the contact and therefore is instru- 
mental in determining the film thickness. Ihe inlet surface and fluid . 
temperatures increase steadily with speed because the film thickness does 
not change appreciably with speed and consequently results in higher 
viscous dissipation at higher speeds. 


27 


B, Ball Surface Temperatiure Measurements Under Severe Conditions 

In order to understand the failure of elastohydrodynamic films and the 
behavior of the contact temperatures, it is necessary to examine the EHD 
contact under conditions which result in relatively thin films, i.e. films 
of the same thickness as the composite roughness of the bounding surfaces. 
Under such circumstances, the fluid film may be locally discontinuous, 
resulting in ambiguous fluid temperature readings. Therefore, during 
this study, only the ball surface temperature was measured. In addition, 
for convenience, measurements were taken only at the center of the Hertaian 
contact. The temperature at the contact center is near the maximum 
(see Figure 7) and is, therefore, representative of the most severe condi- 
tions in the EHD contact. 

Figure 10 shows a plot of the ball temperature rise at the contact 

center above the bath temperature , as a function of sliding velocity for 

2 

peak Herts pressures ranging from 0.52 to 2.03 GN/m . For a given Herts 
pressure, the data plotted on log- log coordinates falls on one. straight 
line for velocities up to a certain value and then on another line, of 
lower slope, for higher velccitiesi This decreasing dependence of tem- 
perature on velocity is supported, by the arguments of Tbrchard [22] and 
Jaeger [23] . Applying Archard’s analysis to a point contact results in . . 
the temperature difference between the fluid midplane temperature Cassumed 
to be the maximum fluid temperature} and the ball siirf ace tempefatu given 
by 
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For a given load, the fluid temperature rise is proportional to the 

2 

product (TC'V-h). At Hertz pressures ^1*5 GK/m f the film thickness and 

traction coefficients were essentially constant, resulting in a fluid 

midplane temperature rise which should be proportional to velocity. It 

is apparent that the ball surface temperature rise is less dependent on 

velocity (Figure 10) and that this dependence decreases with increasing 

velocity. Considering the fluid inidplane as a heat source at temperature 

the ball surface temperature is a function of this temperature level and 

two characteristic times i the time constant for heat penetration to a 

2 

depth of one contact radius below the ball surface t^ is (a Cj^pj^/2kj^) and 

the characteristic fluid resident time t^ is (a/v) [22] . For a given load 

2 

and constant h, which is essentially true for F ^1.5 GN/m , the ratio 
ti/t 2 is directly proportional to velocity. For a given load, in fact, 
t^ is essentially constant under these conditions, which implies that the 
diffusion rate into the ball is fixed, while the fluid resident time 
varies inversely with velocity. A point on the ball surface is then 
exposed for time t^ to the fluid mid plane temperature level, which 
increases linearly with velocity, but the exposure time decreases with 
velocity. The net effect is a surface temperature with a decreasing slope 
as shown in Figure 10 . The temperature rise shown in Figure 10 is the 
difference between the ball surface and lubricant bath temperatures. 

Figure 10 also contains ball surface temperature data as a function 
of load. The same data, at selected velocities, is replotted in Figure 11' 
as a function of Hertz pressure. The temperature rise is the difference 
between the ball temperature at the contact center and the l\f3ricant bath 
temperature as measured vrith a thermocouple. Figure 11 shows that except 
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for the combined conditions of relatively high velocity and high Hertz 
pressure, the data for a given speed follows a power law model with an 
exponent on Hertz pressure of approximately 2. Equation (4) can be 
rewritten in terms of Hertz pressure, resulting in 

-T - T “ P V(TC*h) (5) 

I b 

for a given fluid. The product (TC*h) depends on both P__ and V/ but the 

dependence of the product on is small. Therefore, holding V constant, 

(Tf-Tj^) should be proportional to Herts pressure to the first power. 

However, using the time ratio t^/t^r described above, the larger dependence 

on P of the ball surface temperature rise shown in Figure 11 can be 
H 

quantitatively rationalized. Holding velocity constant, t^/t^ is propor- 
tional to the Hertz radius a which is directly propoarfcional to P . The 
P axis in Figure 11 is also a t^/t axis for a given V. Since it has 
been argued that the heat source has a temperature rise proportional to 
P„ (equation 5) , the temperature rise of the surface, will also depend on 

II 

the relative values of t and t . As P„ increases, t increases as its 

X 2 H X 

square. This means that the effective heat diffusion into the ball is 
decreasing significantly and the surface temperature will remain relatively 
high. In addition, as P„ increases, t- increases proportionately. This 
results in a longer exposure time of the ball surface to the heat source. 
Therefore, although. rises linearly with. Hertz pressure the ball 

surface temperature will rise significantly faster, because of both 
Increasing fluid residence time and decreasing thermal diffusion. 
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c. Comparison of Results with the Blok-Jaeger-Archard Flash Temperatures 
A useful summary of the original flash temperature formulation of 
Blok [24] and the theoretical calculations of Jaeger [23] has been given 
by Archard [22] along with a graphical procedure for determining this 
temperature. In this section, a comparison of ball surface temperatures 
calculated using this method and those obtained experimentally has been 
made as shown in Figure 12. Agreement is good in most cases. The 
calculated values are contact average temperatures predicted by 
theoretical considerations. The measured va]ues are either the 
temperature at the contact center when only such temperatures were 
recorded, or the mean of the highest and lowest temperatures in the 
Hertzian contact area when a temperature map was obtained, A 
regression analysis of measured temperatures on calculated temperatures 
resulted in a coefficient of determination of 0.85. Table II shows 
these temperatures as functions of Hertz pressure and sliding velocity. 

It should be noted that the calculated temperatureo ware obtained by 
adding the bath temperature, as measured with a thermocouple, to the cal- 
culated temperature rise. The measured values, however, whether the 
centerline or contact average temperatures, are the values which result 
directly from the experimental data reduction. The measured values should 
be somewhat higher than the calculated values since both the lubricant 
and surface temperatures just outside the EHD contact are significantly 
higher than the bath oil temperature due to conduction and inlet heating. 

The measured and calculated temperatures agree reasonably well both 
in magnitude and tnrend with increasing speed and load. Because it con- 
siders conduction as the only heat transfer mechanism, the flash tempera^ 
ture theory should best agree with experimental results at high pressures. 
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In these cases f the ratio of contact length to height is the greatest 

thereby making conduction more significant. This trend is supported by 

2 2 

tlie data for P =1.51 GN/m and 1.05 GN/ra , in which contact average 
H 

measured temperatures have been determined. The average deviation 

2 

between calculated and measured temperatures was 4.8% for the 1.51 GN/m 

2 

data and 17% for the 1.05 GN/m data using a Celsius temperature scale. 

Although the flash temperature concept offers a means of estimating 
the contact surface temperature it is subject to substantial limitations. 
The average value obtained can be significantly in error, but more 
important, the method yields only an average. The peak temperature, which 
may be more important in determining contact behavior, is often much 


higher than the average. 
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D. Surface Roughness Effects 

To study the effect of surface roughness on ball surface temperature 
at the contact center, AISI 52100 chrome steel balls of different surface 
roughnesses were used in the EHD contact simulator. Up to this point, 
the center line average surface roughness of the steel balls has been 
. Oil ym. The roughness characteristics of the balls used in all studies 
were independent of direction on the surface. This random orientation 
is a result of different lapping and polishing techniques. 

Profiles of the sapphire surface and ball surfaces were measured 
using a Bendix Group XV System with a modified specimen holder to accom- 
modate the curvature of the ball surfaces. The surface profiles are 
shown in Figure 13 for the sapphire surface and the three ball surfaces 
investigated. In addition to the .011 ]im surface roughness ball, balls 
of ,076 pm and .381 pm surface roughness are shown. The plots in Figure 
13 do not all have the same vertical magnification, as noted. The 
curvature of the ball surface can also be seen at high magnification 
levels. In the following discussions, the .011, .076 and .381 pm surface 
roughness balls will be referred to as the smooth, medium rough and rough 
surface balls respectively. 

Figures 10 , 14 and 15 show ball surface temperatures at the center 
of the EHD contact as a function of speed and load for the smooth, medium 
rough and rough balls respectively. In addition. Figure 16 shows the 
temperature data replotted as temperature rise Cdifference between ball 
surface and bulk oil temperature) as a function of peak Herts pressure 
and sliding speed. In the experiments with the rough ball, the tempera- 
ture reported is the maximum in the contact, whether or not it occurs 



Sapphire Flat 


i 


. 






X' 


2 . 5p in . 3jj m 

T Y~ 


i 

.01 ran 



Medium Rough Ball (.076 ym c.l.a.) 



.01 mm 


Rough Ball (.38 ym c.l.a.) 



Figure 13. Surface Roughness Profiles 


Bill! Temperature 





Tonperature ( C) 



Figure 15. Maximum Ball Surface Temperature - Rough Ball (0.381 um c.l.a.), 
(A < Ir except as noted). 




AT (C) 


Oil c.i.a. surface roughness (smooth), (A > 1) 

076 c,l.a. surface roughness (medium rough), (A < 2) 
381 c.i.a. surface roughness (i^ough) , (A < 1) 


/ /# 

' // / A . 


// 


/ H 

// r 

f/ ' 

V f 

/ '' 

/ / D 



0.2 m/s 
1 m/s 


Max. Hertz Pressure - P„ (GN/m") 

H 


Figure 16. Ball Surface Temperature Rise versus 
Maximum Hertz Pressure. 




41 


at the contact center. It was necessary to make this procediaral change 
because the contact center could not he located when the rough ball was 
used. This is due to the disappearance of the interference fringe pattern 
because of light scattering at the relatively large surface asperities. 
Therefore f the contact was scanned for the maximum surface temperature in 
each case. 

The temperatures for the smooth and medium rough balls are only 
slightly different throughout the range of operating conditions. The 
surface tsnperatures for the rough ball appear to be significantly higher 
than those of the smoother balls. To put all temperatures on an equivalent 
basis the rough ball data in Figures 15 and 16 should be reduced by 
approximately 10°C. This is the difference between the maximum and center 
surface temperatures for the smooth ball (see Figure 3) . 

The parameter A, which is the ratio of EHD film thickness h to the 
composite surface roughness C/. is a recognized parameter for predicting EHD 
contact performance [25]. For values greater than 2 , no asperity inter- 
actions are expected. At A less than 1, severe asperity interaction is 
anticipated. The range 1 ■> A •> 2 is a transition region. Surface rough- 
ness measurements made on the sapphire and ball surfaces in this labora- 
tory are 0 . 1 . a. values. The values of A obtained using c.l.a. rather 
than r.m.s. roughnesses are not sufficiently different that the transition 
values noted above are significantly altered E25J . Also, in computing 
the parameter A, the film thicl^jness value used is that measured for the 
smooth ball. This is necessitated for two reasons. First of ail the 
interference fringe pattern used to determine film thickness disappears 
as the roughness is increased. Secondly, the meaning of film thickness 
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is rjnbiguous at significant roughness levels, since the local 

thickness changes greatly from point to point. The appropriate 

A Values are shown in Figures 10, 14-16. 

In the case of the smooth ball, A is greater than 2 for Herts 

2 

pressures up to 1.810 GW/m , suggesting no asperity interaction. At 
2 

P = 2-03 GN/m A is in the range 1 to 2. However, a subsequent measure- 
H 

ment of the surface profile (the method is discussed in the next section) 

did not reveal any surface alteration. 

Experiments using the medium rough ball resulted in A < 1 for 
2 

P„ > 1,51 GN/m and 1 <; A 2 at lower pressures. Moderate to severe 
H 

asperity interaction is therefore expected at all speeds and loads 

investigated. A relocation surface profile measurement was obtained and 

2 

is shown in Figure 17. At a Herts pressure of 1.70 GN/m it can be 

seen that peaks have been removed in the wear track. 

For the rough ball, A less than one for all operating conditions. 

2 

A profile measurement for P = 1.24 GN/m and V = 5 m/s is shown in 

11 

Figure 18. The removal of asperities is clearly visible on the trace. 

In addition, the alteration of the surface in the wear track could be 
visually detected. 

The above observations reinforce the importance of the parameter A 
in describing the extent of asperity interaction. In addition, it is now 
apparent that significant asperity interaction will result in Increased 

surface temperatures. 
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E, Relo cation Profilimetry 

Figure 19 shows two photographs of the specimen rotation 
attachment with a relocation stage fitted to a Bendix Group XV surface 
measuring system. Because of the curvature of the ball surface, it was 
necessary to modify the standard Bendix system to allow measurements over 
an adequate arc length at high vertical magnification. This has been 
accomplished by rotating the ball about its axis under a fixed stylus. 

A relocation mechanism is incorporated in the rotation attachment. 

In order to detect any change in the surface profile, it is necessary to 
measure the profile at precisely the same point on the ball surface before 
and after being used in the EHD contact. Four of the possible six degrees 
of freedom are removed by using a pair of V -blocks to support the ball by 
two cylindrical collars cemented to the ball surface on opposite sides of 
the ball. The centerlines of the collars and ball are made as nearly 
coincident as possible. Since one of the collars is used to rotate the 
ball in the EHD contact simulator, the plane containing the wear track 
will be normal to this centerline. The other two degrees of freedom 
are removed by providing stops to limit motion along the axis of rotation 
and the amount of rotation about this axis. 

As can be seen from Figure 17, the method has been success- 
ful in relocating the same area on the ball surface for re-examination. 
However, since it is not practical to obtain profiles of the entire wear 
track region before and after running in tbe EHD contact, an indication 
of no asperity interaction is inconclusive. Asperity interaction may 
have taken place at locations other than those measured. 
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P. High Frequency Temperature Fluctuations 

It has been demonstrated in the previous sections that asperity 
interaction can significantly alter the ball surface temperature level. 

The results given thus far, however, are time-averaged temperatures 
obtained with the infrared detector in the DC mode. In this mode of 
operation, the frequency response is 400 Herts. From the surface profile 
measurements (Figures 17 and 18) it has been determined that only a single 
asperity can occupy the detector's field of view (.036 mm diameter) at 
any one time. However, at 1.0 m/s sliding velocity, the residence time is 
only 36 \is. The DC mode, therefore, cannot respond to a temperature rise 
caused by a single asperity interaction. The available AC mode of opera- 
tion, however, can detect such temp«]ii:ature transients. The liquid nitro- 
gen cooled detector has a response time of B ys. The differences in the 
AC and DC modes of operation are shown in Figure 20. 

An important consequence of operating in the AC mode is that the 
reference signal is absent, instead the instrument will produce a voltage 
difference proportional to the variation in target radiation emitted. 

Through an independent experiment, using an external chopper, it has ( 

been determined that the variation indicated in the AC mode is centered 
on the signal received in the DC mode. From this information a plot of 
time-averaged surface temperature along with the maximum and minimum values 
can be obtained. 

Electrical noise problems were encountered when using the AG mode. 

A satisfactory solution was obtained, however, using a variable 
electrical band pass filter. Although the noise and signal could not be 
entirely separated due to the closeness of their frequencies, a pass 




Figure 20. Block Diagram Showing Operation of Baames RM-2A IR Microscope 
in Both DC and AC Modes. 
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band of 1.0 to 20 KHz proved effective, since the high peaks in the 
signal are within 20 KHz for most of the sliding speeds. 

Figure 21 shows a plot of the ball surface temperature at the con- 
tact center as a function of Hertz pressure. The .076 ym c.l.'a. surface 
roughness ball was used at a sliding speed of 1.0 m/s. The plot shown 
was constructed using the DC mode data at five different Hertz pressures. 

In addition, the AC mode was used at tlie four highest Hertz pressure levels. 

As is shown in Figure 21, the AC data shows no fluctuation ^out the DC 

2 

level at P = 1.05 GN/ra , but shows an increasing amount of fluctuation 
H 

as the pressure is increased. The upper and lower curves represent the 
range of temperatures detected. It is believed that the peak values 
represent individual asperity interactions. Figure 21 also shows the sig- 
nificance of the parameter A in predicting the onset of asperity inter- 
action. Figure 22 is a graph of the same data, along with data from 
higher sliding speeds. In all cases, the center value represents the DC 
(time averaged value) , From this data it appears that the range of tem- 
peratures detected in the AC mode of operation increases significantly • 
as more of the normal load is supported by individual asperities . This 
corresponds to a condition of decreasing A values. 



Ball Temperature, 



Figure 21. Average Value and Range o£ Ball Surface Temperature 
Fluctuation - Medium Rougli Ball (.076 pm c. l.a.) , 1 
Sliding Velocity. 
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